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The evolution equations for tensor perturbations in a generic scalar tensor theory of gravity 
are presented. Exact solution are given for a specific class of theories and Friedmann-Lemaitre- 
Robertson- Walker backgrounds. In these cases it is shown that, although the evolution of tensor 
models depends on the choice of parameters of the theory, no amplification is possible if the gravi- 
tational interaction is attractive. 



In spite of the enormous technological problems connected with their detection, gravitational waves (GW) could 
become a very important future source of data in cosmology. Cosmological GW are produced at very early stages 
of the evolution of the universe and decouple from the cosmic fluid at very early times. Consequently they carry 
information about the conditions of the early universe that are not accessible via the electromagnetic spectrum. One 
of the most important applications of a future detection of primordial GW is the possibility of constraining models 
of inflation [l[ . 

Even if GW are decoupled from the cosmic fluid their presence still influences some features of the observable 
universe. In particular, a GW background will produce a signature that can be found in the anisotropies Q and 
polarization [3J of the Cosmic Microwaves Background (CMB). This, together with the extraordinary improvement 
in the sensitivity of CMB measurements, opens the possibility of obtaining important information about GW in an 
indirect way, even though conventional GW detectors are not yet full operational. 

The study of cosmological GW via the CMB is also very important in investigations of alternative theories of 
gravity. In fact, it is well known 0] that the features of GW in GR are rather peculiar and the detection of any 
deviation from this behavior would be a genuine proof of the break down of standard GR. 

For this reason it is important to develop a theory of gravitational waves (or cquivalently of tensor perturbations) 
for alternative theories of gravity. These problems have been widely studied using standard metric-based approaches 
Q, but as yet a fully covariant and gauge- invariant approach, based on the 1+3 formalism 0, @] has not yet been 
developed. The aim of this paper is to show that this method leads to a clear, mathematically well defined description 
of the evolution of tensor perturbations that, unlike metric based methods is easily generalized to non-standard 
gravity. In order to illustrate the power of this approach, we will apply it to one of the most studied extensions of the 
Einstein theory: Scalar Tensor Gravity. This study extends a previous work on the covariant gauge invariant scalar 
perturbations in this framework (8j. 

The paper is organized as follows. In section II we will briefly review the 1+3 covariant approach in a general 
setting. In section III we derive equations that govern the evolution of tensor perturbations for a general imperfect 
fluid. In section IV we adapt these equations the study of scalar tensor gravity and apply them to a specific example. 
Finally, conclusions are given in section V. 

Unless otherwise specified, natural units [h = c = Ub = SnG = 1) will be used throughout this paper, Latin indices 
run from to 3. The symbol V represents the usual covariant derivative and d corresponds to partial differentiation. 
We use the —,+,+,+ signature and the Riemann tensor is defined by 



I. INTRODUCTION 
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where the W a b d is the usual Christoffel symbol (i.e. symmetric in the lower indices), defined by 
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The Ricci tensor is obtained by contracting the first and the third indices 
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Finally the Hilbcrt-Einstein action in presence of matter is defined by 




(4) 
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II. 1+3 COVARIANT APPROACH TO COSMOLOGY 

The starting point (and the corner stone) of our analysis is the 1+3 covariant approach to cosmology pj. This 
approach consists of deriving a set of first order differential equations and constraints for some suitable, geometrically 
well defined quantities (the 1+3 equations) that are completely equivalent to the Einstein field equations. This has 
the advantage of simplifying the analysis of general spacetimes which can be foliated as a set of three dimensional 
(spacelike) surfaces. In the following we give a very brief introduction to the parts of this formalism used in this 
paper. 

A. Kinematics 

In order to derive the 1+3 equations we have to choose a set of observers i.e. a 4-velocity field u a . This choice 
depends strictly on the theory of gravity that we are treating. In this section we give the set of equations for a general 
velocity field and then specialize them to the case of scalar-tensor gravity in section IV. 

Given the velocity u a , we can define the projection tensor into the tangent 3-spaces orthogonal to the flow vector: 

h a b = 9ab + u a u b =>■ h a b h b c = h a c , h ab u b = (5) 
and the kinematical quantities can be obtained by splitting the covariant derivative of u a into its irreducible parts: 

VfcUo = V&Uq - A a U b , Vb1t Q = \Qh ab + C7 ab + UJ ab , (6) 

where V a is the spatially totally projected covariant derivative operator orthogonal to u a , A a = u a is the acceleration 
(a b u b = 0), is the expansion parameter, a ab the shear (a ab — (Tt a b)i aa a = cr ab u b — 0) and u> ab is the vorticity 

(u! ab = W[ ab ], 0JabU b = 0) Hi]. 

In the u a frame, the Weyl or conformal curvature tensor C abcc i can be split into its electric (E ab ) and Magnetic 
(H ab ) respectively: 

E ab = C acbd u c u d =► E\ = , E ab = E [ab] , E ab u b = , (7) 

Hab = 5 Vade C & be U € => H a a = , H a b = H(ab) > H a b U = . (8) 

In what follows we will use orthogonal projections of vectors and orthogonally projected symmetric trace-free part of 
tensors. They are defined as follows: 

v {a) = h a b v b , T< ah > = [ h^ch^d- \h ab h cd ]T cd . (9) 

Angle brackets may also be used to denote orthogonal projections of covariant time derivatives along u a : 

v^=h a b v b , T< Qb > = [ h^ a c h b \~ \h ab h cd }f cd . (10) 

B. Thermodynamics 

The velocity field u a and the projection tensor h ab allow one to split a generic stress-energy tensor T ab as follows: 

Tab = tiu a u b + ph ab + q a u b + q b u a + Tr a b , (11) 

where /i is the energy density and p is the pressure of the fluid, q a represents the energy flux, ir a b is the anisotropic 
pressure and 

q a U a — , ir a a = , TT a b — T^(ab) , ^ ab ^ = ] . 

In a general fluid the pressure, energy density and entropy are related to each other by an equation of state 
p = p(ji,s). A fluid is considered perfect if q a and Tr a b vanish, and barotropic if the entropy is a constant i.e. the 
equation of state reduces to p — p{jS). 
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C. Evolution Equations 

Writing the Ricci and the Bianchi identities in terms of the 1 + 3 variables defined above, we obtain a set of evolution 
equations [25[ : 

leq : rayQ~V a u a = - \ 9 2 + {u a u a ) - 2 a 2 + 2cu 2 - \ (/i + 3p) + A , (12) 

leq : omdotuj^ - \ r] abc V b ii c = - § 6 u a + a% oj b ; , (13) 

leq : sigdota^ - V (a u b > = - § 9 a ab + ii {a ii b) - <r (a c <r b>c - t> w b > - {E ab - \ 7r ab ) , (14) 

{E {ab) + \ ir (ab) ) - (H) ab + \ V (a g fc> = - | (fi + p) a ab - 9 (E ab + 1 7r ab ) (15) 

+ 3cr {a c (EV c - ±7r b > c ) - u^ a 
+ rf d < a [ 2 u c H% + uj c (E% + 1 A) ] , 

F< ah > + (£) ab - I( j7r ) ab = -ei? o6 + 3cr< a c JI b ) c + | w< a q b > (16) 

- r/ cd < a [ 2 u c E% - \ a b \q d - uj c H % } , 

leq : conslfi + V a q a = - G (fl + p) - 2 (u a q a ) - (a a b ir b a ) , (17) 

leq : cons2q^ + V a p + V 6 7r ab = - § 6 q a - a\ q b - (/x + p) u a - u b 7T ab - 7? abc uj b q c , (18) 
and a set of constraints 

leq : divEV b {E ab + ± ^ ab ) - ± V a /z + § 6 q a - \ o\ q b - 3 co b H ab - V abc [ <7 6d ff d c - § co b g c ] = , (19) 

leq : divHV b H ab + (fi + p) co a + 3 u b (E ab - \ 7r ab ) + rf bc [ \ V b q c + a bd (E d c + \ ix d c ) ] = , (20) 

leq : onuV b a ab - § V a 6 + r] abc [ V fc w c + 2 u b w c ] + q a = , (21) 

V a w" - M") = , (22) 

lhconstrH ab + 2 u< a w 6 > + V <a w b) - {<j) ab = , (23) 

that are completely equivalent to the Einstein equations. It is from these equations that we derive the general evolution 
equations for the tensor perturbations in scalar-tensor gravity. 
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III. PERTURBATIONS EQUATIONS 
A. The Background 



The equations presented in the previous section hold in any spacetime we wish to analyze. However, in what follows 
we will focus on the class of spacetimes that do not differ too much from a Friedmann-Lemaitre-Robertson- Walker 
(FLRW) model. The reason is that current observations suggest that the universe appears to deviate only slightly 
from homogeneity and isotropy. We can define a FLRW spacetimes in terms of the variables above. Homogeneity and 
isotropy imply: 

<j = lu = Q, V„/ = 0, (24) 

where / is any scalar quantity; in particular 

V a ^ = V a P = => , u a = . (25) 
It follows that the governing equations for this background are 

e + ie 2 + i(Ai + 3 P ) = 0, (26) 



i? = 2[-ie 2 + M ] , (27) 

fa + e o +p) = o . (28) 

Now in order to describe small deviations from a FLRW spacetime we simply take all the quantities that are zero in the 
background as being first order, and retain in the equations (??-??) only the terms that are linear in these quantities, 
i.e. we drop all second order terms. This procedure corresponds to the linearization in the 1+3 covariant approach 
and it greatly simplifies the system (??-??). In particular, the scalar vector and tensor parts of the perturbations are 
decoupled, so that we are able to treat them separately. In what follows we will focus only on the tensor perturbations. 



B. Covariant decomposition of the 1+3 variables 



Even if the geometrical nature of the kinematical 1+3 variables are known from their definition, the same cannot 
be said for the thermodynamic variables and the terms that contain derivatives. In the first case, this is because they 
depend strictly on the details of the energy momentum tensor of the fluid we are dealing with and in the second case, 
the index carried by the derivative might, in principle, change the geometrical nature of the quantity. For this reason, 
we have to define a general procedure that allows us to recognize the nature of these new quantities or split them 
in their irreducible parts. Giving a completely general procedure is a highly non-trivial task. Fortunately, however, 
for vectors and trace-free symmetric tensor, which describe the key quantities in the covariant approach, this task is 
relatively easy and can be summarized as follows. 

In the case of a vector field the decomposition procedure is based on the Helmholtz's theorem [lOj. This theorem 
says that a continuous vector field V a that vanishes (together with its first derivative) on the boundaries of a manifold 
A4 can be decomposed, modulus a constant, into an irrotational field V a and a solenoidal field V a [26| : 

V a = % + V a = r, abc V b V c + V a v;, (29) 

where 

V a V a = , v abc V b V c = . (30) 

With a little more work we can generalize this decomposition to any symmetric tensor that vanish on the boundary 
of M. [llj obtaining: 

W ab = W ab + W ab + W* ab = W ab + V a W b + V a V b W* , (31) 

where 

V a W ab = 0, {,W)ab = 0, {W*)ab = (32) 
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and both of these decompositions are unique. Now we can define scalars, vectors or tensors as quantities that transform 
like scalars, solenoidal vectors or symmetric tensors, or are obtained from them using the h ab or V a operators. 

In what follows we will consider only purely tensor perturbations which are characterized by the fact that vector 
and scalar contributions to the 1+3 variables and equations vanish, i.e. 

/, W a f, V a V b / = 0, V / (scalar), (33) 

and 

V a , V a V b = 0, V V a (solenoidal vector) . (34) 

C. Tensor perturbation equations 

Focusing only on tensor perturbations, i.e. applying (|33[) and (|34[) . the evolution equations (??)-(??) reduce to 

2 1 

O-afc + - 6 (Tab + E a b - ab = , (35) 

H a b + H ab e + (E) ab -^{ir) ab = , (36) 

Eab +E ab Q- {H) ab + I (p + /Lt) <J ab + ~© TT ab + \-k a b = , (37) 
2 o 2 

together with the conditions 

V b H ab = 0, V b E ab = Q, Hab = {<j)ab- (38) 
Taking the time derivative of the above equations we obtain 

Oab - V 2 CT + ^ 6 Gab + Q O 2 + i (1 - 9 W) f^j G a b = - ab + ^ © ^ab^j , (39) 

H ab - V 2 H ab + 7 - 6 H ab + \ (© 2 - 3 w fx) H ab = (ir) ab + 1 6 i 7T) ab , (40) 

E ab - V 2 Eab+ 7 -QEab + l (<d 2 -3wfi)E ab -lefi (1 + w) (l + 3c 2 ) a ab 
3 3 6 



1 1 - 5 1 



(41) 



where c 2 s is the sound speed of the fluid and we have used the Raychaudhuri equation (??), the energy conservation 
equation ([28]) and the commutator identity 

iX) ab = lX) ab +~(Jt)e. (42) 

These equations generalize the tensor perturbation equations for an imperfect fluid that were derived in [l2j (note 
the different signature!). 

Once the form of the anisotropic pressure has been determined the (|39H4ip can be solved to give the evolution of 
tensor perturbations. As already noted in [l3| the presence of a term that contains the shear in the (|4"Tj) makes this 
equation effectively third order, so that it is not possible to write down a closed wave equation for E ab . If Tr ab — 0, it 
is easy to show that for consistency, the solution for this field must also satisfy a wave equation because the shear is 
a solution of a wave equation and equation (|35[) holds. This will also be the case here because in our case ir ab oc a ab 
and so the anisotropic pressure will also satisfy a wave equation. 
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D. Harmonic analysis 

Following standard harmonic analysis, equations (|39[) and (|40| may be reduced to ordinary differential equations. It 
is standard [Tij to use trace-free symmetric tensor eigenfunctions of the spatial the Laplace-Beltrami operator defined 
by: 

k 2 

V 2 Q afe = — ^Q ab , (43) 

where k = 2ira/\ is the wavenumber and Q a b = 0. Developing a ab and H ab in terms of the Q a b, (|3"5]) and (14"0)) reduce 
to 

+ 1 e *<« + Q e 2 + 1 (i - 9 w) „ - J) ^ = + 1 e tt« , (44) 
ff w + 1 e + \ ( e 2 - 3 w 11 ^) hW = + 1 e c ^) {k) , (45) 

3 3 \ ar ) 3 

and equation (|37|) reads 

£ ;Cfe) = _^W_2 ec7 W + I 7r (fe)_ (46) 



IV. TENSOR PERTURBATION IN SCALAR TENSOR GRAVITY 



A. Scalar Tensor gravity 

In Scalar-Tensor theories of gravity the inclusion of Mach's principle implies the introduction of a scalar field <p 
non- minimally coupled to the geometry [l6| . In particular, this field is associated with the gravitational "constant" G 
and leads to it varying in time. Like many of the extensions of GR, scalar-tensor theories appear in many fundamental 
schemes and have been proposed as models for dark energy because their cosmology naturally leads to the phenomenon 
of cosmic acceleration, which is the characteristic footprint of dark energy [17], [H| ■ 

The most general action for Scalar- Tensor Theories of gravity is given by (conventions as in Wald [lj|) 

A = j dx A J—g [\F($)R - ±V Q 0V> - V{4>) + C m ] , (47) 

where V(<j>) is a general (effective) potential expressing the self interaction of the scalar field and C m represents the 
matter contribution. 

Varying the action with respect to the metric gives the gravitational field equations: 
F{4>)G ab = F{cp) (R ab - § g ab R) = 

= T2 + V a 0V 6 - g ab (±V c </>V c </> + V{^)) + (V 6 V a - 5afc V c V c ) F{<j>) , (48) 

while the variation with respect to the field <f> gives the curved spacetime version of the Klein - Gordon equation 

V a V°0 + iF'(0)i? - V'(<t>) = , (49) 

where the prime indicates a derivative with respect to <j). Both these equations reduce to the standard equations for 
GR and a minimally coupled scalar field when F(<p) = 1. 
Equation (|4"8")l can be recast in the form: 

r < _ rpTOT _ 1 ab , T (eff) f rr.\ 
^ab-^ab — ~pU^j ' ab > l 0U 7 

?d as 

T a e b ff) = T ab = y^) t V ^ V ^ - 9ab (±V c 0V c + V (</>)) + V b V a F((f>) - g ab V c V c F(0)] , 

and an effective gravitational coupling G e // = F(0) _1 appears. Provided that (j) >a ^ 0, equation (j49|) also follows 
from the conservation equations 

V b Tf fc = . (51) 



where T^f^ is defined as 
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B. Perturbation equations in Scalar- Tensor gravity 



At this point we are ready to derive the equations for the evolution of tensor perturbations in scalar-tensor gravity. 
The key point we will use here is that, as we have seen, the field equations of these models can be recast in the 
form (|50p . where together with matter, an effective fluid with energy momentum tensor T^ b appears. Following the 
prescription given in section II, it is then relatively easy to generalize (|39M41[) to this case. 

A crucial point in this procedure is the choice of frame: since (1501) describes a two fluids system, it is not immediately 
obvious which choice of u a is the more convenient. It turns out that if we assume standard matter to be a perfect 
fluid, the best choice of frame is one comoving with the effective fluid. In this frame, standard matter is in general 
imperfect, i.e. it acquires a frame-induced heat flux and anisotropic pressure. However, it can be shown poT ]. that for a 
perfect fluid only the frame-induced heat flux is first order, i.e. standard matter does not produce any frame-induced 
first order tensor contributions. Consequently, in the effective fluid frame we avoid any frame-induced contributions 
to both the effective fluid thermodynamic quantities and standard matter. We therefore obtain the simplest possible 
form for the total anisotropic pressure appearing in the R.H.S. of equations (|39Mlj) . 

In order to define the effective fluid frame, we assume that the momentum density X7 a 4> is timelike [H. Il5l. l2~lj : 



V a ^V Q < , 

in an open region A of spacetime and we choose the 4- velocity to be 

V Q 



l, v = <t> = « a v a = (-v a ^v 



(52) 



(53) 



This implies 



V a = O, 



(54) 



and 



u a b = —h a c h b d \/[d 



= 



(55) 



i.e. our foliation selects vorticity free spacelike hypersurfaces in which <fi = const. 
Using the choice (f53"l) , the tensor T™ T can be decomposed as follows [8] : 



Htot — t L m + 



F{4>) 



S^j 2 + V{<t>) - QF{4>) 



(56) 



PTOT = Pr, 



F(<t>) 



±V 2 - V(cf>) + (F(0) + |6F(0) 



(57) 



TOT 
"a 



9« 



F(<t>) Ua ' 



(58) 



TOT 
T ab 



'ab 



F(d>) 



°~ab ■ 



(59) 



As anticipated, 7r^ does not affect our tensor perturbations equations, so that the total anisotropic pressure in the 
frame u a is proportional to the shear. This allows us to rewrite (|39H4ip as 



j-(fc) 



e 



F 



3 F F 



e 2 



(1 -9w) {!■ 



F 2 



F F 
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I pa g p3 5 pp 1 p 1 p p k 2\ 
~ 6V + 4F-2FF + 6 e F + 2F-^J CT(,=0 - (60) 

These equations describe the evolution of tensor perturbations for a general scalar tensor theory of gravity and reduce 
to the ones presented in [l^, EH in the GR limit (F — * 1). 

In the next subsection, we will apply these equations to a specific example of coupling, potential and background. 



C. The case F(0) = -£<j> 2 , V((p) = W in vacuum 



As an example, let us consider the class of scalar tensor theories given by 

F(4>) = -£(r)0 2 , V{cj ) ) = \^\ 
where £ and p depend on a single free parameter r, via 



£(r) 



(2r + 3) 2 



12(r+l)(r + 2) 



p(r) 



3(1 + r) 
3 + 2r 



(61) 



(62) 



Using the Noether symmetry approach [23j it can be shown that this class of theories admit the following exact 
solution in a vacuum: 



a(t) = at n ; <p(t) = pt m ; K = 0, 



with 



n (r) = 



2r 2 + 9r + 6 
r(r + 3) 



m(r) = — 



2r 2 + 9r + 9 
r(r + 3) 



(63) 



(64) 



This represents a specific case of a general exact analytic solution for the background equations [8, 22, 23|. Depending 
on the value of r, the (|63p can represent power law inflation behaviour (n(r) > 1) or a Friedmann-like phase (0 < 
n(r) < 1) [27j and its character is essentially related to the choice of the form of the potential (via the choice of r). 
The parameter a is free while (3 is linked to A and r through 



A = 



(6 + r)(3 + 2r) 2 
8r(3 + r) 2 (2 + 3r + r 2 ) 



It turns out to be simpler to leave j3 free and derive A, but of course the real free parameter is this last one. Let us 
write down the time dependence of some important quantities. The expansion parameter is 



e(t) = 3 



n(r) 



The effective energy density is given by 



3(6 + 9r + 2r 2 ) 2 1 



a (3 + r) 2 



which is a non negative function of r. Finally the effective barotropic index is 

2 2 + r 

w = — 

3 6 + r(9 + 2r)' 

which is constant. Substituting in the equations (|60|) and (|60l) we obtain 



, (fc) 3(4 + r) (l + 2r) (fc) 
r (3 + r) t 



3 (1 + r) (4 + r) (9 + 2r) 



r(3 



fp.n(r) 



= 



(65) 
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m 24 + 5r(9 + 2r) ^ 
r(3 + r)t 



2(s(2s + 9) +3)(s(5s + 24) + 18) k 



s 2 (s + 3) 2 t 2 



= o 



(66) 



The two equations above give the behaviour of the shear and the gravitomagnetic field. The gravitoelectric field can 
be obtained from the solution of the equations above and 



E (k) 



3 + 9r + 2r 2 



r (fc) 



r (3 + r) i 

In the long wavelength limit (A; — 0), the above equations admit the solution: 



(67) 

*--r , (68) 

(69) 

which are plotted for r = —1 in Figure [U Note that the exponents in the above solution depend on r so that there 
are values of r for which the above solution present growing modes i.e. the gravitational waves are amplified on 
superhorizon scale. Since also the gravitational constant depends on the parameter r, we are interested in values of 
this parameter for which we have growing modes and an attractive gravitational interaction. Unfortunately, it can be 
easily shown that such modes are incompatible with attractive gravitational interaction i.e. the values of r for which 
the exponents in the solutions above are positive imply a negative effective gravitational constant. This means that 
in this class of theories and background only cosmological models driven by a repulsive gravitational constant induce 
an amplification of gravitational waves. 

In the general case, the solutions for (|65[) and (|66p involve Bessel functions of first and second kind J and Y. We 
obtain: 

„W = 



H (k) 




(70) 



(71) 



(72) 



The form of these solutions are similar to the one found in GR and this was expected from the structure of the 
equations. A plot of these solutions for r — — 1 is given in Figure O As in the long wavelength limit also here 
growing modes are possible in principle but not compatible with attractive gravitational interaction so that also in 
short wavelengths only negative G e ft cosmological models induce an amplification of gravitational waves. 



V. CONCLUSIONS 



In this paper we have derived the evolution of the tensor perturbations in a general scalar-tensor theory of gravity. 
These equations have been obtained using the covariant gauge invariant approach to perturbation theory. The key 
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FIG. 1: Plot of the solutions of (|67|) for a scalar tensor theory of gravity in vacuum with F(4>) = — £0 2 , V((f>) = the 
background (I63|) . k = 10 and r = —1. The dashed represents <j^ k \ the dot-dashed curve represents and the solid curve 

represents E^ k ' 



15 r 



10 - 



5 - 




FIG. 2: Plot of the solutions of (|70|l for a scalar tensor theory of gravity in vacuum with F((j>) — ~-(,4> 2 , V{<f>) = \(j) p , the 
background (|63[) . fc = 10 and r = —1. The solid curve represents a (h \ the dashed curve represents and the dot-dashed 

curve represents 



element of this approach is the fact that the Scalar-Tensor field equations can be recast in a form similar to the GR 
field equations, where the energy-momentum is described standard matter plus an effective fluid connected to the 
non-minimally coupled scalar field appears. In this way, equations for the evolution of the tensorial perturbations 
can be derived in a straightforward way. 

To illustrate the method, we analyzed the the evolution of tensor perturbations of a class of theories and backgrounds 
obtained by the Noether symmetry approach [23|. Our analysis reveals that, although the details of the evolution of 
the gravitational waves depend on the specific values of the parameters, no growing modes are compatible with an 
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attractive gravitational interaction. 
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